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Abstract 

The recent result of the higgs search at the LHC experiment has lead to 
more attention to the supersymmetric standard models with heavy sfermions. 
Among them, the models with the almost pure wino being the lightest super- 
symmetric particle (LSP) have been widely discussed due to their success in 
providing a consistent dark matter candidate. The notable phenomenological 
feature of the wino LSP is the degeneracy with its charged SU(2) l partner (the 
charged wino) in mass. The tiny mass splitting makes the charged wino long- 
lived, which allows us to detect the wino production at the LHC experiment 
by searching for the disappearing charged tracks inside the detectors. Since 
the reach of the experiment is sensitive to the mass splitting, it is mandatory 
to estimate it very precisely. We therefore perform a full calculation of the 
mass splitting at two-loop level, and find that the splitting is reduced by a few 
MeV compared to the one-loop calculation. This reduction leads to about a 
30 % longer lifetime of the charged wino, with which the current constraint on 
the wino mass by the ATLAS experiment is improved by about 10%. 



1 Introduction 



The supersymmetric standard model (SSM) is one of the most attractive candidates 
for physics beyond the standard model (SM). Both the discovery of the higgs bo- 
son [U [2] and the null-observation of supersymmetry (SUSY) signals at the LHC 
experiment have given us some hints for SUSY model buildings. In the minimal 
SSM (the MSSM), for example, the higgs boson mass is predicted to be smaller 
than the Z boson mass at the tree level. The observed mass of the higgs boson 
at around 126 GeV, therefore, indicates that huge radiative corrections to the higgs 
self-coupling from the SUSY breaking effects are required [3] — [S] - 

One of the simplest scenarios leading to such huge corrections is putting the 
masses of the sfermions (especially of the squarks) at the scale of 0(10-100) TeV [3] 
[7J . It should be noted that although the squarks are far beyond the accessible range 
of the LHC experiment, in such cases, this class of the scenarios does not necessarily 
mean that all the SUSY particles are as heavy as 0(10-100) TeV. For example, if 
we suppose that the SUSY breaking field is charged under some (gauge) symmetries, 
gauginos cannot acquire their masses through the linear term of the SUSY breaking 
field in the gauge kinetic functions of the MSSM. In this case, the leading contri- 
butions to the gaugino masses come from the anomaly mediated contribution [HJ |9], 
which are one-loop suppressed compared to the squark masses. The gaugino masses 
are therefore predicted to be (9(0.1-1) TeV which are accessible at the LHC ex- 
periment. This class of the high-scale SUSY scenarios has recently attracted more 
attention, and phenomenological and cosmological aspects of the scenarios have been 
studied extensively pTO] — [23] J 1 

One of the most prominent features of the anomaly mediated gaugino spectrum 
is that the lightest supersymmetric particle (LSP) is predicted to be the almost pure 
neutral wino, which is highly degenerate with the charged wino in mass. This is 
because the mass splitting between these two particles is forbidden at the tree-level 
due to the approximate custodial symmetry. The mass splitting is dominated by the 
radiatively generated contributions which are estimated to be 160-170 MeV at the 
one- loop level [25] [27]. 

1 In this paper, we base our discussion on the pure gravity mediation scenarios |10| . where the 
Higgsinos and the heavier Higgs bosons in the MSSM are predicted to be as heavy as the sfermions. 
Our formulas obtained in the following analysis are applicable to any heavy sfermion scenarios as 
long as the Higgsinos are much heavier than the gauginos as in the scenarios such as PeV-Scalc 
Supersymmetry [Mj and Spread Supersymmetry jTTJ . 
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Due to the degeneracy between the charged and the neutral wino masses, the 
charged wino decays mainly into the neutral wino and a soft pion with the decay 
length of cr = (9(1-10) cm, which allows the wino production to be detected by 
looking for a disappearing charged track inside the detectors at the LHC experi- 
ment [2B]- This signal is characteristic for the high-scale SUSY scenarios with the 
anomaly-mediated gaugino mass spectrum. The expected number of the charged 
track is proportional to exp[— L/cr] with L being the distance between a detector 
and a collision point, and hence, the reach of the experiment is very sensitive to cr. 
Therefore, a precise calculation of the decay length is mandatory. 

In this article, we calculate the mass splitting between the charged and the neutral 
winos at the two-loop level. In Ref. [29], the splitting of the winos has been evaluated 
at the two-loop level in the heavy wino limit by calculating non-decoupling contribu- 
tions. The result, however, cannot be directly applied to the wino mass in the range 
of 0(100) GeV where the LHC experiment is searching for the winos. We there- 
fore perform a full two-loop calculation of the splitting including the non-decoupling 
effects. In the next section (section[2]), we calculate the mass splitting at the two- 
loop level. We will see that the contributions from the SM particles are dominant, 
while those from the SUSY particles are negligible. In sectional we discuss the decay 
length of the charged wino and compared with the recent experimental results by the 
ATLAS collaboration [28]. As a result, we find that the decay length of the charged 
wino becomes 10-30 % longer than that obtained at the one-loop calculation, with 
which the current constraint on the wino mass by the ATLAS experiment turns out 
to be severer than the LEP2 constraints [30] [33]. Section H] is devoted to summary 
of our discussion. 



2 The mass splitting 

As mentioned above, the neutral wino (x°) and its charged SU(2)^ partner (the 
charged wino, x^) are almost degenerated in mass at the tree-level due to the ap- 
proximate custodial symmetry. The dominant mass splitting, 5m = m^± — m^o, is 
generated by radiative corrections, which pick up the breaking of the custodial sym- 
metry as pointed out in [27]. In this section, we calculate the radiative corrections 
at the two-loop level. 
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2.1 SM contributions 

When the sfermions, Higgsinos, and the heavier higgs bosons are in the range of 
(9(10-100) TeV and decouple from the low energy physics below the TeV scale, the 
neutral and the charged winos only couple to the SM particles through the SU(2)^ 
gauge interaction. In such cases, the radiative correction to the mass splitting from 
the SM sector can be calculated by using the effective Lagrangian, 

-9 (xW + /i-c.) + g X ~ {cwZ + s w 4) X", (1) 

where £sm is the SM Lagrangian and M 2 is the invariant mass of the winos. The 
notation for the SM gauge fields is understood, and SU(2) L gauge coupling is denoted 
by g, while cw(sw) = cos6'vi/(sin6 l iy) with 9w being the weak mixing angle. 

The mass splitting between the charged and the neutral winos is caused by the 
custodial symmetry breaking by U(l)y gauge and Yukawa interactions. It should be 
noted that the breaking of the custodial symmetry is highly suppressed at the tree- 
level in the wino-SM system. In fact, at the tree-level, the breaking of the custodial 
symmetry is mediated through the Higgsino mixing. As a result, the tree- level mass 
splitting is highly suppressed by the Higgsino mass, fi, which is given by, 

m^(sin2/3) 2 tan 2 #ty 14 keV / 300 GeV \ / 100 TeV\ 2 
Imixing ~ (Mi - M 2 )/i 2 ~ T^J \M 1 - M 2 ) V ) ' { ) 

Here, mw denotes the mass of the jy-boson, (3 the higgs mixing angle of the MSSM, 
and Mi the mass of the bino, 2 ! As we will see below, the above tree- level mass 
splitting is sub-dominant compared to the radiatively generated mass splitting]^ 

2.1.1 The pole mass 

The pole mass of a spin half particle can be extracted from the 1PI effective two-point 
function, 

r 2 = f - M + Z K (p 2 )p + S M (p 2 ) , (3) 



2 The mass splitting in Eq. ([2]) is valid for Mi — 3> mz- 

3 In the Split Supersymmetry models [33] where the Higgsino can be as light as the gauginos, 
the tree-level contribution to the mass splitting is not necessarily negligible. 
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with p being the four momentum of the particle and Mq the tree-level mass. Thus, 
for given self-energy functions, E# and E^, the pole mass is iteratively given by, 



M 



pole 



Re 



M n 



Sm(M p 2 o1 J 



(4) 



l + ^(M p 2 ole ) 

In a perturbative analysis, we expand the above pole mass as a power series of 
coupling constants. At the two-loop level, the above iterative expression of the pole 
mass is reduced to 



M, 
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Mi — E 
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Here, E^ M and E^ M are the self-energy functions at the one- and two-loop levels, 



respectively, while the dotted functions, E 
respect to p 2 . 



(i) 



r, denote the derivatives of E^ M with 



2.1.2 Renormalization scheme and input parameters 

We take the input parameters to the above effective Lagrangian; 

a, rh w , rhz, M 2 , m u m hl and Q, 

where the hatted variables denote the MS variables, and Q is the renormalization 
scale. All the quark and lepton masses except for the top quark mass are neglected 
in our analysis. 

To relate the above listed input parameters (the MS variables) to the experi- 
mental observables, we have to take finite renormalization effects into account. In 
the following analysis, we extract the input parameters by using the renormalized 
relations at the one-loop level, 



a 



SM 
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i+m(Q 2 )/Q' 
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m\ - Il zz (m 2 z ) 



M 2 -M 2 ESJ ) (m|o)-E^(m| ) , 



(6) 
(7) 
(8) 
(9) 



where all the self-energies (H xx and E^ M ) used in the above equations are given 
in the appendix [B] Here, dg^ denotes the QED fine structure constant in the MS 
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W+ w° w + w + w° w ± 



Fi gure 1: One-loop diagrams contributing to the functions k * n 

scheme in the SM at the Z-boson mass scale, mw,z the physical W and Z boson 
masses, m^o the physical neutral wino mass. It should be noted that the one-loop 
relations are precise enough for the two-loop estimation of the wino mass splitting, 
since the leading mass splitting starts at the one-loop level. 

The top quark and the higgs boson appear only at the two-loop calculation of 
the mass splitting. Thus, the MS variables m t and m^ may be replaced with their 
physical masses m t and rrih at this level of precision. As for the top quark mass, 
however, we use the MS top mass at the one-loop level for m$ As we will see, the 
Q dependence of the mass splitting at the two-loop level comes mainly from those 
of the top mass m t . We set, on the other hand, m/, = m/, since the running of the 
higgs mass does not cause significant effects on the splitting. 

Once we obtain the input parameters, d, m w , and m z from Eqs. ©-flE]), we 
can calculate g, g' using tree level relations. In deriving the one-loop relations in 
Eqs. (jS}-©, we also obtain the counter-terms to subtract ultra-violet (UV) diver- 

(2) 

gences. These counter terms play important roles to calculate H K M , as will be 
discussed later. 

2.1.3 The mass splitting at one-loop level 

The one-loop result of the mass splitting between neural and charged winos is well 
known [25] [27] and used in many literatures. The loop diagrams of the winos and 
gauge bosons shown in Fig.[T] lead to the functions S^ M . With the use of the 
formula in Eq. (E]) and the self-energies given in the appendix[B] the mass 

splitting 5m = m^± — m^o at the one-loop level is given by 

= (^ 2 M 2 /87r 2 )[/(myM 2 2 )-c 2 w /(m|/M 2 2 )], (10) 

4 The finite renormalization effect connecting between rht (MS mass) and mt (pole mass) is the 
same as those in the SM, because the scalar top quarks are heavy and decoupled. 
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where the function f(z) is defined as f(z) = J Q dx(l + x) log[l + z(l — x)/x 2 }. In the 
heavy wino limit, M 2 ^> rnz,w, the mass splitting is reduced to 

5m ~ —(m w - c 2 w m z ) , (11) 
which is about 160-170 MeV. 

2.1.4 The mass splitting at two-loop level (strategy) 

The two-loop self-energies, S^ M (M|), are obtained from the two-loop diagrams 
(Fig. [2]) and from the diagrams including counter-terms which cancel the one-loop 
UV divergences (Fig.[3D • In our actual analysis, we first calculated the two-loop 1PI 
amplitudes using FeynArts [35] and FeynCalc [36], which were reduced to a set of 
basis integrals by TARCER [37]. We finally evaluated the integrals numerically using 
TSIL [38]. For the diagrams including counter-terms, we used the ones given in 
the appendix IB. 31 As a nontrivial cross check, we have confirmed that all the UV 
divergences are properly canceled. 

We also have to care about infra-red (IR) singularities. For the charged wino, 
the amplitude in Fig.[2]-(i) in which a photon is circulating in the outer loop and the 
one in Fig.[3J-(b) with the photon loop behave as 

and hence, they are IR divergent. In addition, the derivatives, ^k]m\ p 2 =m^ are a ^ so 
IR divergent due to the diagram including a photon propagator. We have checked 
that all the IR divergences are canceled with each other when we evaluate the pole 
mass in Eq. fl5]). See the appendix O for more discussions on the cancellation of the 
IR divergences. 

2.1.5 The mass splitting at two-loop level (result) 

Now, let us show the resultant mass splitting at the two-loop level. In the followings, 
we take 6t^{m z ) = 127 '.944 ± 0.014, m w = 80.385 ± 0.015 GeV, m z = 91.1876 ± 
0.0021 GeV [39], m t (m t ) = 163.3 ± 2.7 GeV @0] and m h = 125.5 ± 0.7 GeV as the 
SM input parameters^ 

5 The ATLAS collaboration reports the Higgs boson mass as 125.2±0.3±0.6GeV [41], and CMS 
collaboration reports 125.8 ±0.4 ±0.4 GeV [12]. Here, we take naive average over these results, and 
combine statistical and systematic error. 
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(d) (e) (f) 




(c) (d) 

(2) 

Fi gure 3: Diagrams including counter-terms which contribute to the function in 
Eq. |5]). The counter-terms are determined to renormalize one-loop divergences. 
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(a) m neutra | ino =100 GeV 
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(b) m neutralino = 1000 GeV 
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Figure 4: The renormalization scale dependence of 8m. The green lines show 5m 
at one-loop level in Eq. 07J)) . and the red lines is 5m at the two-loop level which is 
evaluated by Eq. (EJ) in MS scheme. We take m^o =100 GeV (a) and 1000 GeV (b). 
Here, we take m t {m t ) =163.3 GeV and m^ =125.5 GeV. 



S 



In Fig.lH we show the renormalization scale dependence of 5m, which is the 
dominant source of the theoretical uncertainty of the mass splitting at the two-loop 
level. In the figure, the two-loop result is numerically evaluated by Eq. (JSJ) in the MS 
scheme, while the one-loop result is evaluated by Eq. (fTOj) in terms of the on-shell W 
and Z boson masses, i.e., 

5m Uoop = [f(m 2 w /Ml) - c^/(m|/M 2 2 )] , (13) 

where we defined c^, = m^/m z and + = 1. The Q dependence of the one- 
loop result in Eq. (I13p comes from the running of the gauge coupling constant, while 
it comes from all the MS parameters appearing in by Eq. (EJ). The figure shows that 
the Q dependence becomes weaker at the two-loop level as expected, since the mass 
splitting should not depend on Q at full-order. In our analysis, we found that the 
dominant source of the Q dependence of the two-loop result is the running of the 
top quark mass. 

Since the uncertainty of the mass splitting due to the choice of Q is expected to be 
compensated by the three-loop contributions including the QCD and the top- Yukawa 
interactions, the Q dependence of the two-loop result gives us a rough estimation of 
the uncertainty of the mass splitting from the higher-loop effects. In the followings, 
we fix the renormalization scale as Q = m t at which the mass splitting is least 
sensitive to the renormalization scale Q, and estimate the uncertainty of the mass 
splitting due to the choice of Q by, 



d5m 
dlogQ 



(14) 



-m t 



In addition to the above uncertainty, there are expected to be other uncertain- 
ties from the higher-loop corrections which are not encapsulated in the choice of 
the renormalization scale. At the three-loop level, for example, the dominant non- 
decoupling contribution to the mass splitting is expected to be proportional to m t 
and the QCD coupling. Although the numerical factors of those corrections cannot 
be determined unless explicitly calculated, we give naive estimations to those higher 
loop corrections by, 

^s-ioop^m = (|^) 2 [^j 7im t ~ 0.033 MeV. (15) 

Here, we have multiplied a factor of tt which is expected to accompany the non- 
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Type of error 


Estimate of the error 


Impact on 8m 


asM(mz) 


experimental uncertainty in «sm(^z) 


±0.018 MeV 


mw 


experimental uncertainty in mw 


±0.019 MeV 


m z 


experimental uncertainty in mz 


±0.001 MeV 


rh t 


experimental uncertainty in m t 


±0.081 MeV 


m h 


experimental uncertainty in m^ 


±0.002 MeV 


Experiment 


Total combined in quadrature 


±0.085 MeV 


choice for Q 


QCD and top Yukawa at one- loop by Eq. (fl4l) 


±(0.3 - 0.4) MeV 


three-loop 


naive estimation by Eq. (fT5l) 


±0.033 MeV 


Theory 


Total combined in quadrature 


±(0.3 - 0.4) MeV 


Total 


Total combined in quadrature 


±(0.31 -0.41) MeV 



Table 1: Experimental and theoretical errors in the evaluation on 5m at two-loop 
level. 



decoupling effects at M 2 ^> m$ 

The experimental errors of the input parameters also lead to uncertainties of the 
mass splitting. As we summarize in Tab.fU however, the effects of the experimental 
errors are relatively small compared to the theoretical errors. As a result, we find 
that the uncertainty on 5m is dominated by the three loop logarithmic corrections, 
i.e. the renormalization scale dependence. 

In Fig.0 we show the mass splitting between the neutral and the charged winos 
as a function of the neutral wino mass|j The figure shows that the two-loop contri- 
butions reduce the mass splitting by a few MeV compared to the central value of 
the one-loop result. We also show the theoretical and experimental uncertainties as 
green/red bands. As a result, we find that the uncertainties are significantly reduced 
by the two-loop analysis. For the sake of readers, we give a fitting function of the 



We have confirmed that the naive estimation of the two-loop contribution. 



a 2 N 2 



A 2 -ioo V Sm = (-±J nm t ~ 3.9 MeV , (16) 

gives a fair estimation of our two-loop numerical results. 

7 We have also confirmed that the limit mw,z *C m^, our result reproduces the one in Ref. [29] 
at this level of precision in the heavy wino limit, M% S> mz- 
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central value of the two- loop result for Q = m t , 
8m , /, 



-413.315 + 305.383 flog ) - 60.8831 flog 

V 1 GeV/ V 



"X" 



1 MeV ' V to 1 GeV/ ' V to 1 GeV/ 

+5.41948 f log - 0.181509 flog r ^ v j . (17) 

for the central values of the SM input parameters. Deviation of the above fitting 
function from our two-loop result is smaller than 0.02 % for the wino mass being 
100-4000 GeV. 



2.2 SUSY contributions 

Before closing this section, let us evaluate the contributions to the mass splitting from 
the diagrams including the heavy SUSY particles in the tens to hundreds TeV range. 
Since the winos couple to the other gauginos (bino and gluinos) only through the 
exchanges of those heavy particles, all the SUSY contributions to the wino masses can 
be expressed by the higher- dimensional operators suppressed by the heavy masses. 
At a first glance, a five dimensional operator, 

C 5 = je abc x a x\HW c H) , (18) 

with H being the light higgs boson and A = 0(10-100) TeV the scale of the heavy 
SUSY particles seems to break the custodial symmetry and contribute to the mass 
splitting. Here, the superscripts a, b and c denote the indices of the adjoint represen- 
tation of SU(2)^. This operator, however, vanishes because of the Majorana nature 
of the winos, x a x b = X b X a - Another dimension five operator 

£ 5 = ±x a (H^T a H)b , (19) 

with b being the bino, on the other hand, contributes to the mass splitting of 
0(v 4 /A 2 Mi) through the neutralino mass matrix. Incidentally, the tree-level mass 
splitting due to the Higgsino mixing in Eq. ([2]) can be regarded as one of the contri- 
butions of this type with A ~ fi. As a result, we find that the contributions from 
the dimension five operators are negligibly small as we have seen in Eq. ([2]). 

The next lowest dimensional operator which contributes to the mass splitting is 
the dimension seven operators 

£r = ^-Ax a X b ){H ] r a H){H^r b H) , (20) 
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170 




Figure 5: The wino mass splitting 5m as a function of m^o. The dark green 
band shows 8m at the one-loop level which is evaluated by Eq. / [TO]) with uncertainty 
induced by Q dependence, and the red band shows 5m at two-loop which is evaluated 
by Eq. ([3J) in MS scheme. The light green band shows the uncertainty for one-loop 
result evaluated by Eq. ( fl6j) . The uncertainties for the two-loop result induced by the 
SM input parameters and the non-logarithmic corrections are negligible (see Tab.\Tj). 
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where M denote the insertion of the gaugino massE For A = (9(10-100) TeV, the 
contribution from this operator to the mass splitting is again negligibly small. 

3 The charged wino decay 

As we have seen in the previous section, the charged and the neutral winos are 
highly degenerated in mass. Therefore, the decay width of the charged wino is highly 
suppressed by the phase space integral, and hence, the charged wino is long-lived 
and has the decay length about cr = (9(1-10) cm. With such a rather long decay 
length, it is possible to detect the charged wino production at the LHC experiment 
by looking for disappearing tracks. In this section, we estimate the lifetime of the 
charged wino and compare with the constraint from the disappearing track search 
by the ATLAS collaboration [2%] . 

With the small mass splitting 5m ~ 160 MeV, the charged wino dominantly 
decays into a neutral wino and a soft charged pion. At the leading order, the decay 
width of the charged wino can be expressed in terms of the decay width of the 
charged pion, 

T( X ± -»• X°vr ± ) = r(vr ± fuj x 1 - -JL i _ _£ , 21 

m^m^ \ om z J \ J 

where m n and denote the masses of the charged pion and the muon, respectively]^ 
The branching ratio of the sub-leading leptonic decay mode into a pair of the electron 
and the neutrino [H] , 

7C 2 

rXx^fWe) * Y5^ 6m5 - (22) 

We consider the above two decay modes. 

In Fig.lHl we show the lifetime of charged wino as a function of the charged wino 
mass, m^±. The meaning of the green and red bands are the same with the ones 



8 This operator can be obtained from, for instance, a dimension eight operator 
(qLX a H)^ (qLX a / A 4 which is generated by integrating out the squarks (especially stops) at the 
tree-level. By integrating the quark-loop and inserting the gaugino mass, we obtain the dimension 
seven operator in Eq. (|20l) . 

9 At the next-to-leading order, Eq. (f2T|) receives radiative corrections from the QED and the 
electroweak interactions which are expected to be around (a/ir) log(m^/m v ) ~ 2 %. In this paper, 
we neglect these corrections to the total decay width and leave the detailed analysis of the decay 
width for future study |43j . 
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Figure 6: The lifetime of charged wino evaluated by using 5m at the one-loop (green 
band) and two-loop (red band). We neglected the next-to- leading order corrections 
to the lifetime of the charged wino estimated in terms of the pion decay rate, which is 
expected to be a few percent correction. We also neglected the leptonic decay mode 
of the charged wino which branching ratio is about 1 %. The black chain line is the 
upper limit on the lifetime for a given chargino mass by the ATLAS collaboration 
at 95 % CL (y^i = 7 TeV, C = 4.7 fb _1 ) [2EJ - 
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in Fig.[5j The region above the black chain line is excluded by the disappearing 
charged track search by the ATLAS collaboration at 95 % CL The figure shows 
that the lifetime is enhanced by about 30 % for the wino mass around 100 GeV due 
to the two-loop contributions. Furthermore, the figure also shows that the precise 
estimation of the mass splitting at the two-loop level improves the constraint on the 
charged wino mass by about 10%. 

4 Summary 

We have calculated the mass splitting of the charged and the neutral wino in the 
almost pure wino LSP scenario at the two-loop level. Such a scenario attracts more 
attention after the discovery of the Higgs-like boson at the LHC experiment. As a 
result, we found that the lifetime of the charged wino becomes about 30% longer 
due to the two-loop contributions. Furthermore, we also found that the precise 
determination of the mass splitting greatly improve the constraint on the mass of 
the wino mass obtained by the disappearing track search at the ATLAS experiment 
by about 10%. 
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A Loop functions 



The one-loop functions are presented in this appendix, which are used to calculate 
radiative corrections to the self-energies of SM particles and winos; 

„ , 9 9 9x . f 1 , , (1 — x)m 2 + xmk — x(l — x)p 2 — ie , . 

B (p 2 ,m 2 ,m 2 ) = A-j^dxlogl >—k L_J >£ , ( 23 ) 

BiOAmLml) = -j + J^dxxlogl ^ ^ , (24) 

A Z" 1 21 {I - x)m\ + xm 2 - x{\ - x)p 2 - ie 

B 2l {p 2 ,m{,m 2 2 ) = --/ dxx 2 log^ i * ^ , (25) 

3 Jo Q 

where A = 2/(4 — d)— 7^+log(47r). With M being much larger than m, the functions 

B (M 2 ,M 2 ,m 2 ) and B X {M 2 , M 2 , m 2 ) are evaluated as follows; 

B (M 2 ,M 2 ,m 2 ) = A-log^ + 2-^ + (^log^) , (26) 

rw„.2 ,,2 2n A 1, M 2 3 vrm ^ / m 2 , M 2 \ 

fll (i#»,««,m') = -- + - 1 „ g _-- + - + (_l„ g _] . (27) 

By using B (p 2 , m 2 , m 2 ,), -E>i(p 2 , m 2 , m 2 ,), and B 2 i(p 2 , m 2 , m 2 ,), we define HV(p 2 , m 2 , ml 
and B 22 (p 2 ,m 2 ,ml) by, 



n v (p 2 ,m 2 ,m 2 ) = -p 2 [B 1 { P 2 ,ml,ml) + B 2l (p 2 ,m\,ml)) , (28) 

9 2 1 



o 2 1 

£ 22 (p 2 ,m 2 ,m 2 ) = V(^i+^2i)-^S --(m 2 -m 2 )(S + 2 J B 1 ) . (29) 



B Radiative corrections at one-loop 

Here, all radiative corrections to the 1PI self-energies of the gauge bosons and the 
winos at one-loop level are presented. The counter-terms (in MS scheme) to eliminate 
the one-loop UV divergences are also shown. These self-energies as well as the 
counter-terms are used in the calculation of the mass splitting at two-loop level. 

B.l Gauge boson self-energies 

In terms of the 1PI amplitude IT(g 2 ), the full propagator (2-point function) of the 
gauge boson in the Feynman gauge is given by, (—ig^u)/[<l 2 — m>y + Il(g 2 )]. In this 
subsection, we present the contributions to the amplitude II (q 2 ) from both the SM 
particles and the winos, which are divided into three parts; 

u VlV2 = n^g + + ng } V2 + P 2 5 ZviV2 + , (30) 
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where V1V2 = 77, 7Z, ZZ, and WW. The first term in the right-hand side is the 
contributions from the quarks and the leptons, the second term is those from the 
gauge-higgs sector of the SM, and the third term is from the neutral and charged 
winos. The fourth and fifth terms show the counter terms given in appendix IB. 31 

B.l.l Contributions from winos 



c 2 



n $V) = ^2n v (p 2 , M|, M|) , (31) 

n$V) = -^-n y (p 2 ,M|,M|), (32) 

n zi(P 2 ) = ^C-n^M^Mf), (33) 

n w(P 2 ) = J^n y (p 2 ,Mf,Mf). (34) 

With the use of the above amplitudes, the finite renormalization effect, 11$ (p 2 ), in 
Eq. (jSJ is given by the combination, fl$(p 2 ) = H$(p 2 ) — p 2 (e 2 /127r 2 )A. 

B.1.2 Contributions from quarks and leptons 



<TV) = -E^2^^/ n v(P 2 >/>/) , (36) 

(37) 



/ 

n 2 N C 

U ZZ IP J - 2^ 27 r2r 2 
j Z7r C W 

n 2 N C 



T 2 \ m 2 
-^- + Z 2 ) U v (p 2 ,m},m 2 ) - -^B (p 2 , m 2 , m 2 ) 



,2 



U v (p 2 ,ml,m 2 d ) - ^Bi{p 2 ,rh 2 u ,m 2 d ) 



(38) 



where Qf is the electric charge of the fermion /, while Tf takes the value 1/2 and 
— 1/2 for up-type fermions (u, c, t quarks and neutrinos) and down-type fermions 
(d, s, b quarks and charged leptons), respectively. The coefficient Zf is given by the 
equation Zf — (Tf/2 — Q/s 2 ^), while iVj? = 3 for the quarks and Nj? = 1 for the 
leptons. The summation ^2f ,^ should be over left-handed quarks and leptons. The 
finite mass effect of the down- type fermions on the amplitude TlfywiP 2 ) * s neglected. 
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B.1.3 Contributions from the gauge-higgs sector 



<§(p 2 ) 



3e^ 
~4vr 2 



6/2-2 - 2 \ , P 2 

B22(P ,m w ,m w ) + — 



^ 2 R (r? ™ 2 rh 2 \ 



+ 



eg ,~ 2 ins / 2 -2 - 2 \ i egcwP 2 
-ip—iS>c w - 1)B 22 {V ,m w ,m w ) + 

r C\y Z47T 

Q o g ( 2d wP 2 + m 2 v )B {p 2 ,m 2 v ,m 2 v ) > 

I^r *»(p >%.%) " 

/ *2 2,o-2 -2\n / 2 -2 -2 \ 

~Tp~2~( WP +2m w -m z )B (p ,m w ,m w ) 

t 

g 2 

-TT^i^wP 2 + 3 ™w - rh 2 z) B o(p 2 , rhw> rh%) 

107T 

<? 2 ~ 



(39) 



(40) 



[B22(p 2 ,m 2 z,m 2 h ) ~ m 2 z B (p 2 ,m 2 z , rh 2 h )\ , (41) 



(42) 



B.2 Wino self-energies 

With the use of the 1PI amplitudes Sj<-(g 2 ) and £m(<? 2 ), the full propagators (2- 
point functions) of the winos are given by, + Y* K (q 2 )}(ji — M 2 + T, M (q 2 )}. In this 
subsection, we explicitly present the amplitudes for both neutral and charged winos 
at the one-loop level. For the neutral wino, the amplitudes are given by 

^2 f n 



E 
£ 



(i) 



K,0 
(1) 



g 

16tt 2 
g 2 M 2 r 



4 J B 1 (p 2 ,M 2 2 ,m 2 ^) + 2 
8B {p 2 ,M 2 ,m 



J M,0 16n 2 

On the other hand, the two amplitudes for the charged wino are given by 

s 2 w B 1 (p 2 Jd 2 ,0)+c 2 w B 1 (p 2 ,M 2 ,m 2 z ) + B 1 (p 2 ,M 2 ,m 2 w ) + l] + S z ~ 
s 2 w B (p 2 , Ml 0) + c 2 w B (p 2 , M 2 ,m 2 z ) + B a (p 2 , M 2 ,m 2 w ) - l] - 5 M , 



v (i) 9 



g 2 M 2 



J M,± 



4tt 2 



where explicit forms of the counter-terms, 5z- and 5m-, are given in IB. 31 



(43) 
(44) 

(45) 
(46) 
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B.3 Counter terms 

Finally, we give the counter-terms in the MS scheme in the framework of the SM 
plus the winos. These are used in the calculations of the self-energies mentioned 
above and of the mass splitting at the two-loop level as shown in Fig.[3l 

B.3.1 Gauge boson self-energies 



5 Zzz 
$Zww 



32 

16^ 2 V 9 9 



A 



eg 



I6ir 2 cw 

\2 



4 32s 2 w 
3 + n I "9 



9 



16vr 2 



16vH 



c 2 
w 



-N, 



9 



(t 






3 )\ 



A 



% + ~FT S W I ^ + 



11 H.2 5 -4 

Y + y Sw ~ 3 Sw 



where JV„ is the number of the generation, namely N„ — 3 for the SM. 



^m 2 

m 7 Z 



6 m% z 



eg 



IQtt 2 cw 



16ir 2 c 2 w 



(2-2s 2 w )m 2 z A , 
3 

-m 2 + (-1 + Qs 2 w - 4s^)m| 



A , 



16vr 2 



-mi + {-l + 2s 2 w )m z 



A, 



(47) 
(48) 
A , (49) 
(50) 

(51) 
(52) 

(53) 



where, we have neglected the masses of all the SM fermions except the top quark. 
B.3. 2 Wino self-energies 



6* 



8tt 2 



g 2 M 2 

' 2tt 2 



A 



(54) 
(55) 



B.3. 3 Gauge interaction of the wino 

The neutral and charged winos have the SU(2)x, gauge interaction which is described 
by the term, £ int = ie a b c (g + 5^w)x a ^^ rb X C i an d the counter term is given by 

,3 



5 



xxW 



9 
Air 2 



A . 



(56) 
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C Cancellation of IR divergences 



In Eq. (ED with M being M 2 , the IR divergences appear in [(e£ ) +M E^ ; )(2M E^ 
2M 2 E^)] p 2 =M 2 and (— + MqE^- 1 ). The first term is the products of the one- 
loop contributions. The one- loop amplitude (E^ + M S^) is explicitly written 
by 



W) = E&V)+M>e£V) = | 



ri 4 fc (ie 2 )7^(^+^ + M ) 7M 
(2vr) 4 [fc 2 -m2][(A;+p) 2 -M 2 ; 



+ 



(57) 



where we have introduced a photon mass m 7 to control the IR divergences. The 
ellipses stand for the contributions from the loop diagrams of the W and Z bosons, 
which are nothing to do with the IR divergences. The derivative of the one-loop 
amplitude Fi^q 2 ) with respect to p 2 gives the IR divergent contribution, 



_d 

dp 



2 iW 



(2vr) 4 [k 2 -m 2 ][{k+pf 



M 2 } 2 



+ ••• , 



(58) 



where the ellipses represent the terms which do not cause the IR divergences, namely 
the IR-safe terms. 



The second contribution, F 2 l 
written by 



-1)(E 



(2) 
M 



(2)> 



is, on the other hand, 



F 2L (M 2 ) = (*e 2 ) J 



d 4 k -ftf + M )[E«(A; 2 )^ + E$(* 2 )](tf + Mohp 



[(k - p) 2 - m 2 ][k 2 - M 2 f 



+ 



.(59) 



) 2 =M„ 2 



The numerator of the integrand in above equation is possible to be simplified as 

r2M Q [M^\k 2 ) + E$(fc 2 )](j5 + Moh» + ^{k 2 - M 2 ) 
= 2M [MoE^ (M 2 ) + E$ (M 2 )]2M o7 ^ + 0(k 2 - M 2 ) . (60) 

As a result, the IR-divergent part of the two-loop contribution can be reduced to, 



(2M ) \M V%\M$) + E%(M$)\ x ^F 1L (p 2 



d 



+ 



(61) 



Therefore, we find that the IR-divergences cancel with each other, and hence, the 
pole mass is an IR-safe quantity. 
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